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The possibility of self-consistent determination of instanton liquid parameters is discussed to- 
gether with the definition of optimal pseudo-particle configurations and comparing the various 
pseudo-particle ensembles. The weakening of repulsive interactions between pseudo-particles is 
Qh' argued and estimated. 
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The problem of finding the most effective pseudo-particle profile for instanton liquid (IL) model 
^ . of the QCD vacuum pQ has already been formulated in the first papers treating the pseudo-particle 
^! superposition as the quasi-classical configuration saturating the generating functional [2] of the fol- 
Q-T lowing form 

j§ Z= f D[A] e~ s ^ , (1) 
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where S(A) is the Yang-Mills action. Although the solution proposed in Ref. [2] was quite acceptable 
\ phenomenologically the consequent more accurate analysis discovered several imperfect conclusions 
t^j- ' putting into doubt the assertion about the instanton ensemble getting stabilization and some addi- 
tional mechanism should be introduced to fix such an ensemble [3]. In this note we revisit the task 
formulated in Ref. [2] within the self-consistent approach proposed in our previous paper [I] . We are 
not speculating on the detailed mechanism of stabilizing and are based on one crucial assumption 
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t**- , which is the existence of non-zero gluon condensate in the QCD vacuum. This idea is not very orig- 
^ 1 inal but turns out far reaching in the context of our approach. The particular form and properties 

of this condensate will be discussed in the following paper. 
^ Thus, as the configuration saturating the generating functional ([1]) we take the following super- 

^ , position 

A a „(x) = BZ(x)+J2A«(x ni ) , (2) 

i=l 

here A ^ stands for the (anti-)instanton field in the singular gauge 

Al{x- 1 ) = -u a %^% f{y), y = x-z, (3) 
9 V 

7i = (pi,Zi,Ui) denotes all the parameters describing the z-th (anti-)instanton, in particular, its size 
p, colour orientation u>, center position z and as usual g is the coupling constant of gauge field. 
The function f(y) introduces the pseudo-particle profile and will be fixed by resolving the suitable 
variational problem. For example, for the conventional singular instanton it looks like 



f(y) = — 7j2 ■ (4) 
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In analogy with this form we consider the function / depending on y 2 or, more precisely, on the 
v 2 

variable x = *jp at some characteristic mean pseudo-particle size p. Dealing with the anti-instanton 

one should make the substitution of the 't Hooft symbol fj —> tj. It is seen from ([2]) we 'singled out' 
one pseudo-particle of ensemble and introduced the special symbol B for its field which actually has 
the same form as Eq. (jHJ). 

The strength tensor of this 'external' field and the field of every separate pseudo-particle A can 
be written as 

G% = G%(B) + G%(A) + G%(A, B) , (5) 
where two first terms are given by the standard definition of field strength 

G%(A) = d,A a u -duAl + g r hc A\Al , (6) 

with the entirely antisymmetric tensor f abc . In particular, for the singular instanton of Eq. (jHJ) it 
takes the form 
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where /' means the derivative over y 2 . The third term of Eq. (JHJ) presents the 'mixed' component of 
field strength and is 

G^(A, B) = g f*°(B>AZ - B b v AD = g f a V ~ (Sj f)du« ~ B b u ffa) ^ /. (8) 

It was shown in Ref. [1] that in quasi-classical regime which is of particular interest for appli- 
cations, the generating functional (JTJ) could be essentially simplified if reformulated in terms of the 
field Bj± averaged over ensemble A. Performing the cluster decomposition [5] of stochastic exponent 
inEq. (pQ) 

where (Si) = ((Si)), (S1S2) = (Si)(S2) + ((S1S2)), ■ ■ ■ (the first cumulant is simply defined by 
averaging the action) the higher terms of effective action for the 'external' field in IL could be 
presented as 



and the mass m is defined by the IL parameters developing for the standard singular pseudo-particles 
(HI) the following form (see, also below) 

m 2 = 9n 2 np 2 jj^, (11) 

with n = N/V where N is the total number of pseudoparticles in the volume V and iV c is the number 
of colours. The small magnitude of characteristic IL parameter (packing fraction) np 4 allows us at 
decomposing to keep the contributions of one pseudo-particle term (~ n) only. 

The effective action in Eq. (1101) implies a functional integration in which the vacuum stochastic 
fields are not destroyed by the external field. Then there is no reason to develop the detailed 
description of the field B driven by the symmetries of initial gauge invariant Lagrangian for the 
Yang-Mills fields. In practice it could be understood as an argument to do use the averaged action 
dealing with the field B. It means the colourless binary (and similar even) configurations only of 
field B survive in the effective action. In other words the decomposition B ~ B^ + • ■ ■ is used 
(in what follows we are not maintaining the index for the field B). Obviously, if there is any need 
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of more detailed description including, for example, information on the fluctuations of field B one 
should operate with the correlation functions of higher order and the corresponding chain of the 
Bogolyubov equations. 

The selfconsistent description of pseudo-particle ensemble may not be developed based on Eq. 
ffTUj) only because in such a form the pseudo-particles of zero size p = are most advantageous. 
In Ref. [I] the version of variational principle was proposed which makes it possible to determine 
the selfconsistent solution in long wave-length approximation for the pseudo-particle ensemble (anti- 
instantons in the singular gauge with standard profile PJ) and external field. Here it adapts to the 
saturating configuration also and its more optimal (than standard) profile is defined, as suggested 
in Ref. [2], taking into account the IL parameter change while the pseudo-particle field is present. 

The contribution of saturating configuration into the generating functional is evaluated as (see 
[2J for the denotions) 

r 1 r N 

Z^Y= D[B]— ft d 7i e -s(*r> . (12) 

i=l 

The following terms should be taken into consideration 

N N 

S(B, T ) = - E ^ d (pi) + P U ™t + E u Lt(B) + S(B) , (13) 
i=i i=i 

(the details of deducing this expression can be found in [1]). Here we remind only that to obtain it 
one should average over the pseudo-particle parameters and to hold the highest contributions only 
at summing up the pseudo-particles. If the saturating configurations are the instantons in singular 
gauge with the standard profile (jlj) the first term describing the one instanton contributions takes 
the form of distribution function over (anti-)instanton sizes 

,b b-5n2N, 



where 



j3 = -Mn(Ap), 



d(p) = C Nc A b p b - b (3 2 ^, (14) 



b = T Nc ~ l Nf ' (15) 



4.66 exp(-1.68iV c 



7r 2 (iV c -l)!(iV c -2)! ' 

If one considers the profile of Eq. the change of one pseudo-particle action which has the form 

S, = 3 r*£/3 \{y 2 ff + f{l-ff] , (16) 

should be absorbed while calculating. Here (3 = 8n 2 /g 2 is the characteristic action of single pseudo- 
particle (jlj) which is defined at the scale of average pseudo-particle size (3 = (3(p) where (3{p) = 
— lnCN c —bhi(Ap). The coefficient b enters the corresponding equations (in particular the distribution 

function ( TT4T) ) always with the additional factor s = It means that in all the formula containing 

the one instanton contribution the following substitution 

b^bs. (17) 

should be done. The penultimate term of Eq. ffl3|) accumulates the partial pseudo-particle contribu- 
tions coming from the 'mixed' component of the strength tensor (jSj) and describing the interaction 
of pseudo-particle ensemble with the detached one, i.e. 



uUb) = I d A x 



G^(A h B) G^jA, B) 



3 



The other terms at the characteristic IL parameters are small as it was shown in Ref. jl]. The 
average value of 'mixed' component is given by the following formula 



1 8 M 1 2 f 2 

KM, B) G%(A, B)U = VW ^T I K Sj , B 2 = - 2 J - 2 , (11 

v i \ c i y y 



here / is defined by the integrated profile function of pseudo-particle 

/VcyJjB 7T 2 O 2 f°° y 2 

dy -^-r- f , I = —— / dx f 2 , x = — . 
4 Jo p z 

In particular, for the standard form of pseudo-particle we have 

The corresponding constant (see [1]) Co = ^-f~ n*-i should De changed for the modified one 

POO 

C = ACo , A = / dxf 2 , 
Jo 

in all terms describing the interaction of IL with detached pseudo-particle if the profile function 
/ is arbitrary. Eq. (fl8l) demonstrates that we are formally dealing with non-zero value of gluon 
condensate which is given by the correlation function 



(Al(x; i)Al(y\ i)) u , = j 2 -^-^ F {^—^ ) . (19) 



For the pseudo-particle of standard form the function F(A) equals to 



7T 2 A 2 + 2 



F(A) = — ^ VA 2 + 4 In 



VA 2 + 4(A 2 + 1) + A 3 + 3A 



VA 2 + 4 - A 

2 ^ r V , a2\ , 2 A 2 



A 2 

with the asymptotic behaviours 



__ (A 2 + l) 2 ln(1 + A2) + ^ A , ^ | A 



(20) 



lim F(A) -> tt 2 - ^ A 2 + vr 2 A 2 In |A| , lim F(A) - £ . 



The presence of this condensate ( 1191 which leads, in particular, to the mass definition as in (fill just 
signifies the assumption mentioned at the beginning this note. 

The second term of ffl3l) describes the repulsive interaction between the pseudo-particles of en- 
semble 



P Uint =J2j d 



4 x I ^-yX^A?) ^JAi^j} 



hi \ ' 7i>7j 

and actually presents the same contribution as U ex t but being integrated with the field B of every 
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individual pseudo-particle as (3 Ui nt = J d 4 x B 2 . It results in the change of coupling constant 
£o = 27 4 2 Ji-i describing the pseudo-particle interaction (see [2]) for new form 

e = a 2 e , 

(similar to the change of constant (). And eventually the last term of Eq. (fl3l) presents simply the 
Yang-Mills action of the B field 



S(B) = J d A x 



It is worthwhile to notice that the topological charge of the configuration (Tj0) is retained to be equal 
to 

1 r G a G a r°° 1 

N = I d " X = " 6 Jo dX /7(1 " /} = 1 ' °^ = 2 °^ ' 

here e^ UQ/ 3 is an entirely antisymmetric tensor, e 12 34 = 1. 

The generating functional (fl2l) might be estimated with the approximating functional (see [2]) as 

Y>Y 1 exp(-(S-Si)) , (21) 

where 

r 1 r N 

Yi = J D[B] — / J] e-*™- s W , S 1 (B, 1 ) = -^ln^Pi) , 

i=l 

and p(p) is an effective one particle distribution function defined by solving the variational problem. 
In our particular situation the average value of difference of the actions is given as follows 

1 1 r N 

{S ~ Sl) = YwJ II d lt [pU mt + U ext ( 1 ,B)-J2^d( Pt ) + ^^)}e^ hlM = 

N r a(p) 8 N 2 1 r 

= — dp p(p) In-— + - 772 — d ll d l2 Uintijx,^) p(pi)p(p 2 ) + 
p J d(p) 2 V Pq J 



+ f*» j dp p 2 C B' 2 

J V J Un 



with pq = J dp p(p). In this note we estimate the functionals in the long wave length (adiabatic) 
approximation, i.e. consider the IL elements to be equilibrated by the external fixed field B. After- 
wards, with finding the optimal IL parameters out we receive the effective action for the external 
field in the selfconsistent form. Eq. (1221) is taken just in such a form in order to underline the inte- 
gration is executed over the IL elements and the parameters describing their states are the functions 
of external field (i.e. could finally be the functions of a coordinate x). The physical meaning of such 
a functional is quite transparent and implies that each separate IL element develops its characteristic 
screening of the attached field. 

Now calculating the variation of action difference (S — Si) over p(p) we obtain 

= c d( P ) e -(^v+cB 2 )P 2 s 

where C is an arbitrary constant and its value is fixed by requiring the coincidence of the distribution 
function when the external field is switched off (B = 0) with vacuum distribution function then 

p( p ) = C N J 2N <A bs p bs - 5 e-Wef+CB*)? (23) 
With defining the average size as 

— [dp P 2 Ml) 

P = , 

Po 

we come to the practical interrelation between the IL density and average size of pseudo-particles 

(n{3e7+(B 2 )7^v, (24) 

where v = ^ s 7^ - . Apparently, the size distribution of pseudo-particles can be presented by the 
well-known form as 

H(p) = C N J 2Nc A bs p bs ~ 5 e p 2 . (25) 
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Figure 1: The energy E(a) when the profile function includes a screening effect (1291) with the pa- 
rameter A (s — 1) only taken into consideration (lower curve) and with both parameters used (upper 
curve) (see the text). 



Eqs. ( 1221 and (|23|) allow us to get the estimate of generating functional (12"TT) in the following form 

Y > f D[B] e" 5(B) e~ E , (26) 



E 



n 



d x n I In — -r 

A 4 



2 2 



C Nc f3 



2N C 



Now taking into account Eq. (1241) and fixing a field B, parameters s and A the maximum of functional 
(126p over the IL parameters can be calculated by solving the corresponding transcendental equation 



(^p = 0) numerically. Here it is a worthwhile place to notice the presence of new factor in the 

denominator of k what is caused by the Gaussian form of the corresponding integral over p 
squared and, hence, the integration element requires the introduction of 2p dp. In Ref. [2] this factor 
was missed. However, this fact has not generated a serious consequence because any application of 
these results is actually related to the choice of suitable quantity of the parameter A entering the 
observables (the pion decay constant, for example). It means we should make the proper choice of 
basic scale. Besides, we should also keep in mind the approximate character of IL model. Further we 
give the results for both versions to demonstrate the dependence of final results on the renormalized 
constant Cn c - 

Searching the optimal configuration / we take the effective action in the form of nonlinear func- 
tional as 



S, 



eff 



d x 



E[B\ 



(27) 



in which the IL state is described by solutions p[B, s, A], n[B, s, A]. In practice the following differ- 
ential equation should be resolved 



d 2 f 



1 df , f{l 



2„,2 



d 2 y 



V dy 2 



+ 



/)(l-2/) 1 dE 
6/3o df 



y 1 



(28) 



at fixed initial magnitude of /(x ) putting up the derivative in the initial point /'(x ) in such a way to 
have the solution going to zero when x is going to infinity Parameter /3 is introduced to fix a priori 
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Figure 2: The IL density as the function of x = y 2 /p 2 . Three dashed curves correspond to the 
different profile functions. The lowest dashed line corresponds to the standard form (j3J). The top 
dashed line corresponds to the profile function with the screening factor (I2TJ1) and one parameter A 
(s = 1) included and the middle line presents the same function but with two parameters included. 
The solid line presents the selfconsistent solution of variational problem. 



unknown value of coupling constant in the pseudo-particle definition If the profile function has 

been fixed the configuration should be found in the form in which the starting values of parameters 

s, A and (3q coincide (within the given precision) with the parameters obtained from the solution 

/. Nowadays this approach looks the most optimal one among other existing possibilities not only 

because of the computational arguments but in view of the poor current level of understanding the 

interrelation between perturbative and non-perturbative contributions while calculating the effective 

Lagrangian. In fact, it was mentioned in Ref. [2] that in more general (realistic) formulation of 

this problem Eq. ( 1281) should include the term responsible for the change of 'quantum' constant 

Cn c with the function / changing. In principle, it could imply that the problem of pseudo-particle 

ensemble stabilization is connected at the fundamental dynamics level with the anticipated smallness 
dC 

of the ^ c contribution and, apparently, should be addressed not so much to the description of 

the interacting pseudo-particles and their interactions with the perturbative fields but rather to 
investigation of the time hierarchy corresponding to the breakdown of quasi-stationary behaviour 
of the vacuum fluctuations which will certainly lead to the changes of suitable effective Lagrangian 

In order to receive the preliminary parameter estimates we consider the simplified model with 
the profile function containing only one additional parameter for describing the screening effect as 
regards 

-aX 2 

— = (29) 

The energy E as the function of the screening parameter a is depicted in Fig. 1. The lowest dashed 
curve shows the behaviour when the changes related to weakening of repulsive interaction are taken 
into account by switching on the parameter A only (at s = 1). The top dashed curve was obtained 
with both parameters switched on. The optimal value of the screening parameter a is determined by 
the minimum point of function E(a). Besides, this figure demonstrates the stability of variational 
procedure of extracting the IL parameters. For the first calculation the values of characteristic 



7 



parameters for corresponding solution were taken as a = 0.06, A = 0.775, s = 1.0067 with the 
following set of the IL parameters pA = 0.3305, n/A 4 = 0.919, (3 = 17.186. These values give for the 
ratio of average pseudo-particle size and average distance between pseudo-particles the quite suitable 
quantity p/R = 0.324. For another calculation we have treated the parameter set characterizing 
the solution as a = 0.02, A = 0.888, s = 1.0015 and for the IL parameters the following values 
pA = 0.315, n/A 4 = 0.829, (3 = 17.67, p/R = 0.3. In order to get more orientation we would like to 
mention that for the ensemble of standard pseudo-particles (a = 0, A = 1, s = 1) the corresponding 
values are pA = 0.301, n/A 4 = 0.769, (3 = 18.103, p/R = 0.282. 
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Figure 3: The average size of IL pseudo-particles as the function of x = y 2 /p 2 . Three dashed curves 
correspond to different profile functions. The lowest curve corresponds to the standard form The 
top dashed curve corresponds to the profile function with the screening factor (I2"9"j) which includes 
one parameter A (s = 1) and the middle line shows the same function with two parameters included. 
The solid curve corresponds to the selfconsistent solution of the variational problem. 

Now we examine the impact of correction introduced in Eq. (j26p when we changed the term h ' 
which has been obtained in Ref. [2]. For the first calculation with the set of solution parameters as 
a = 0.24, A = 0.546, s = 1.029 we have for the IL parameters pA = 0.331, n/A 4 = 1.844, (3 = 17.173 
which lead to the ratio discussed equal to p/R = 0.386. For another calculation we have the following 
results a = 0.05, A = 0.799, s = 1.0053 and pA = 0.291, n/A 4 = 1.356, (3 = 18.483, p/R = 0.314. 
And for the ensemble of standard pseudo-particles (a = 0, A = 1, s = 1) these parameters are 
pA = 0.265, n/A 4 = 1.186, (3 = 19.305, p/R = 0.277. 

The Fig. 2 and Fig. 3 show the behaviours of IL density and average pseudo-particle size as the 
functions of distance x. The dashed lines on both plots correspond to the similar ensembles. The 
lowest curves demonstrate the behaviours for the ensembles of standard pseudo-particles (jl]). The 
top curves present the ensemble of pseudo-particles with the profile function (|29|) at a = 0.06 and 
s = 1. And the middle dashed lines correspond to the profile functions with a = 0.02 and s ~ 1.03. 
Obviously, it may be concluded that including even small change of the second parameter value 
(s ~ 1.03) leads to the noticeable change of ensemble characteristics (for example, the IL density) 
because the highest contribution to the action when the coupling constant becomes the function of 
p is essentially modified. 

Let us make now several comments as to the 'complete' formulation of the problem of analyzing the 
equation (|2"H1) . It was numerically resolved by the Runge-Kutta method. This approach combined 
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with numerical calculation of the derivative -jj- at every point of consequent integration interval 

allows us to avoid the problems which appear when searching the minimum of complicated functional 
in multidimensional space. 

v 2 

The initial data were fixed at the point xo = jp- = 0.1. Since the IL density value at the 

coordinate origin is inessential the initial form of pseudo-particle profile function is taken without 
any deformations as /(x ) = ^ j - . Then at fixed values of the parameters A, s and (3q the coefficient 
c is calculated. It allows to set the slope of trajectory /'(x ) = — c/(l — /)/x at initial point in such 
a form in order to have the solution going to zero at large distances. Afterwards we find out the 
values of parameters A and s requiring the input data to coincide with the output ones within the 
fixed precision. The parameter values which obey the imposed constraints are the following (input 
values) A = 0.69099, s = 1.049, f3 = 16.26 at c = 1.361 and A = 0.691, s = 1.049, f3 = 16.263 
(at the output of variational procedure). The solid line in Fig. 4 shows the obtained profile / as 
v 2 

the function of x = >±%. The differences of profiles are smoothed over if they are presented as the 

functions of y because the large magnitude of the screening coefficient, for example a = 0.06, is 
compensated by enlargening the pseudo-particle size. The dashed lines on this plot show the profile 
functions for the standard form (J3J) (top dashed line), with the screening factor (T2U]) including one 
parameter only a (s = 1) (lowest dashed curve) and two parameters included (middle dashed line). 

1 .0 1 1 1 1 1 1 1 1 1 1 1 1 i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i 




Figure 4: The various profile functions. The top dashed curve corresponds to the standard form (j4l), 
the lowest dashed curve shows the function with the screening factor (I29p including one parameter 
A (s = 1) and the middle line presents the same function with two parameters included. The solid 
line corresponds to the selconsistent solution of variational problem. 



Another calculation (with modified T-function contribution) was based on the slightly different set 
of relevant parameters which are for the input values A = 0.607, s = 1.0515, /?o = 17.04 at c = 1.545 
and A = 0.6066, s = 1.0515, j3o = 17.042 for the output one at the finish of variational procedure. 
The behaviours of IL density and average pseudo-particle size for selfconsistent solution are plotted 
in Fig. 2 and Fig. 3 (solid lines, respectively^. In the Table 1 we present the IL parameters at 

It is interesting to notice that considering IL (ensemble of pseudo-particles in the singular gauge) in the field of 
regular pseudo-particle we obtain the IL density value in the center of regular pseudo-particle which is larger than its 
value at large distances what looks like the anti-screening effect. 
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the large distances from pseudo-particle (the first line) together with the data for the ensemble of 
pseudo-particles with the standard profile function (the second line). The third and fourth lines of 
this Table 1 are devoted to the calculations with the second set of parameters (with factor 2 absent in 
Eq. (|26p ). The fourth line, in particular, presents the calculations for pseudo-particles with standard 
form of profile function. 

Table 1. Parameters of IL. 



pA 


n/A 4 


P 


p/R 


-4 

np 


0.381 


0.743 


16.263 


0.354 


1.582-10~ 2 


0.331 


0.769 


18.103 


0.282 


6.277-10 3 


0.354 


1.245 


17.042 


0.379 


1.955-10~ 2 


0.265 


1.186 


19.305 


0.277 


5.849-10- 3 



It is quite obvious that the utilization of optimal pseudo-particle profile function leads to the larger 
pseudo-particle size but the packing fraction parameter holds, nevertheless, a small quantity which is 
quite suitable for the perturbative expansion. Besides, the results obtained allow us to conclude that 
with tuning A a fully satisfactory agreement our calculations of pseudo-particle size, the ensemble 
diluteness and gluon condensate value with their phenomenological magnitudes extracted from the 
other models are easily reachable. The calculations of several dimensional quantities in our approach 
are also very indicative. The values of the screening mass (|lip . average pseudo-particle size and IL 
density obtained for two values of A (200 MeV and 280 MeV) are shown in Table 2. The sequence of 
line meanings is identical to that in Table 1 as well as the meanings of last four lines which present 
the results of calculations with the second set of parameters (with factor 2 absent in Eq. ( 1261) ). 

Table 2 . Screening mass and IL parameters 



A MeV 


m MeV 


p GeV- 1 


n fm 4 


200. 


381 


1.906 


0.7496 




304 


1.503 


0.7688 


280. 


533 


1.361 


2.88 




426 


1.074 


2.95 


200. 


456 


1.77 


1.245 




333 


1.325 


1.186 


280. 


638 


1.264 


4.78 




466 


0.946 


4.56 



Another interesting feature of this calculation is the weakening of pseudo-particle interaction. This 
effect is driven by the coefficient £ 2 (~ A 2 ). Our estimates for the first set of parameters give 
A = 0.691 and, hence, A 2 ~ 0.48 and for the second set we have (A = 0.607) and A 2 ~ 0.37. Let us 
mention here that the reasonable description of instanton ensemble can be reached in the framework 
of two-component models [6] as well. 

Our calculations enable us to conclude that dealing with IL model (formulated in one-loop ap- 
proach) one is able to reach quite reasonable description of gluon condensate even being constrained 
by the values of average pseudo-particle size and other routine phenomenological parameters. More- 
over, the ensemble of pseudo-particles with standard profile functions turns out to be very practical 
because introducing the other configurations to make the similar estimates is simply unoperable. 
With such an approximation of the vacuum configurations the coefficient of interaction weakening 
develops the magnitude about A 2 ~ 0.3 — 0.5. Including this effect leads to the enlargening of 
pseudo-particle size. It allows us to conclude that nowadays the instantons in the singular gauge is 
the only serious instrument for effective practising. 

The authors are sincerely grateful to A.E. Dorokhov and S.B. Gerasimov for interesting discussions 
and practical remarks. The financial support of the Grants INTAS-04-84-398 and NATO PDD(CP)- 
NUKR980668 is also acknowledged. 
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